LOGARITHMS OF ITERATION MATRICES, AND PROOF OF A 
CONJECTURE BY SHADRIN AND ZVONKINE 
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Abstract. A proof of a conjecture by Shadrin and Zvonkine, relating the en- 
tries of a matrix arising in the study of Hurwitz numbers to a certain sequence 
of rational numbers, is given. The main tools used are iteration matrices of 
formal power series and their (matrix) logarithms. 



This note is devoted to the study of the somewhat mysterious-looking sequence 

II ^ 11 3 1129493 2711 12406 2636317 
' '~2' 2'~3'12'~4'~ir'T'T2"' 6~ ' 15~' 60 '"' ( ' 

of rational numbers. I first encountered this sequence in ongoing joint work with 
van den Dries and van der Hoeven on asymptotic differential algebra [4]. It also 
appears in a conjecture made in a paper by Shadrin and Zvonkine [31] in connection 
with a generating series for Hurwitz numbers (which count the number of ramified 
coverings of the sphere by a surface, depending on certain parameters like the degree 
of the covering and the genus of the surface). I came across [31] by entering the 
numerators and denominators of the first few terms of (S) into Sloane's On-Line 
Encyclopedia of Integer Sequences [1]. (The numerator sequence is A134242, the 
denominator sequence is A134243.) In this note we prove the conjecture from [31]. 
In the course of doing so, we identify a formula for the sequence (S): denoting its 
nth term by c„ (so c\ — 0, C2 = 1, C3 = — | etc.), we have 

r_-n*+i 



lsjfc<n 
Kn 1 <---<n k - 1 <n k = 



Hrww.j n k 

ni\\n 2 \ \rik-i 



Here and below, we denote by {^} the Stirling numbers of the second kind: is 
the number of equivalence relations on a j-element set with i equivalence classes. 
They obey the recurrence relation 

n fi " 11 +J i : 1 i a,j>o) 



(i,j>0). 



i I I i — 1 I i 



with initial conditions 
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For example, we have 

1(311 1 

1 - - < > = 1 3 = — = c 3 

2 (2 J 2 2 



2 V 1 2 J 13J/ 3 [2J 1^3 J 2 V ' 3 2 

1 " * ({& + {£} + ti» + 1 f 1-1(15 + 25+10)- 



I ({^m + + {X}) - = |(3 • 25 + 7 • 10 + 6 • 10)- 

immm U -3-e. io 
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A key concept for our study of (S) is the iteration matrix of a formal power series; 
these matrices are well-known in the iteration theory of analytic functions [20, 21] 
and in combinatorics [11]. The iteration matrix of a power series / S Q[[z]] of 
the form / = z + z 2 g (g € Q[[z]]) is a certain bi-infinitc upper triangular matrix 
with rational entries associated to /. After stating the conjecture of Shadrin and 
Zvonkinc in Section 1 and making some preliminary reductions, we summarize some 
general definitions and basic facts about triangular matrices in Section 2 and intro- 
duce the group of iteration matrices in Section 3. In Section 4 we determine its Lie 
algebra of infinitesimal generators, by slightly generalizing results of Schippers [30]. 
These results tie in with a notion from classical iteration theory: the infinitesimal 
generator of the iteration matrix of a formal power series / as above is uniquely 
determined by another power series itlog(/) G rQp]], introduced by Jabotinsky 
[21] and called the iterative logarithm of / by Ecallc [13]. Some of the properties of 
iterative logarithms are discussed in Section 5, before we return to the proof of the 
conjecture of Shadrin-Zvonkine in Section 7. The exponential generating function 
(egf) of the sequence (c„), that is, the formal power series 

turns out to be nothing else than the iterative logarithm of the power series e z — 1. 

The iterative logarithm itlog(/) of any formal power series / satisfies a certain 
functional equation found by Jabotinsky [20]. In the case of / = e z — 1, this 
equation leads to a convolution formula for Stirling numbers (and another formula 
for the terms of the sequence (c n )): 



(_l)fc+l J- n2 1 j- n3 -| J- Uk 



l<ni<-"<fifc_i<rifc=n 



k [ni J {n 2 J {n k - 

(-l) k (n 2 \ (n 3 \ ( n k 



l^k<n— 1 
Kni<---<nk-i<nk=n— 1 



(C) 



To our knowledge, this formula does not seem to have been noticed before. (For 
instance, it does not appear in Gould's collection of combinatorial identities [17].) 
We give a proof of (C) in Section 7. 

Shadrin and Zvonkine write that the sequence (S) seems to be quite irregular 
[31, p. 224]. This impression can be substantiated as follows. A formal power series 
/ G C[[z]] is said to be differentially algebraic if it satisfies an algebraic differential 
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equation, i.e., an equation 



P(*,/,/',...,/< n >) = 



where P is a non-zero polynomial in n + 2 indetcrminates with constant complex 
coefficients. The coefficient sequence (/„) of every differentially algebraic power 
series / = X)«>o e QIH] i s regular in the sense that it satisfies a certain kind 
of (generally non-linear) recurrence relation [28, pp. 186-194]. A class of differ- 
entially algebraic power series which is of particular importance in combinatorial 
enumeration is the class of D-finite (also called holonomic) power series [32, Chap- 
ter 6] . These are the series whose coefficient sequence satisfies a homogeneous linear 
recurrence relation of finite degree with polynomial coefficients. Equivalently [32, 
Proposition 6.4.3] a formal power series / € C[[z]] is D-finite if and only if / satisfies 
a non-trivial linear differential equation 



(This class includes, e.g., all hypergeometric series.) In Section 7 we will see that 
the egf of (c„) is not differentially algebraic. This is a consequence of a result of 
Boshcrnitzan and Rubel, stated without proof in [10], which characterizes when the 
iterative logarithm of a power series satisfies an ADE; in Section 6 below we give a 
complete proof of this fact. It is also known [8, 25] that the egf of (c„) has radius 
of convergence 0. Indeed, a common generalization of these results holds true: 
the egf of (c„) does not satisfy an algebraic differential equation over the ring of 
convergent power series. The proof of this fact will be given elsewhere [3] . It seems 
likely (though we have not investigated this further) that the ordinary generating 
function (ogf) 



of the sequence (S) is also differentially transcendental. (Note, however, that there 
are examples of sequences of rationals whose egf is differentially transcendental yet 
whose ogf is differentially algebraic; see [26, Proposition 6.3 (i)].) 

Notations and conventions. We let d, m, n, k, possibly with decorations, range over 
N = {0,1,2,...}. All rings below are assumed to have a unit 1. Given a ring R we 
denote by R x the group of units of R. 

Acknowledgements. We thank the anonymous referees whose corrections and sug- 
gestions improved the paper. 



Before we can formulate this conjecture, we need to fix some notation. Let K be 
a commutative ring and let R = K[[to,ti, . . .]] be the ring of powers series in the 
pairwise distinct indeterminates t ,ti, . . . , with coefficients from K. We equip R 
with the m-adic topology, where m is the ideal (to,t\, . . .) of R. In this subsection 
we let i, j range over the set of sequences i = (io,h, ■ ■ ■ ) G N N such that i n = 
for all but finitely many n. For each i we set 



oaf + aif + ■■■ + a n f {n) = (a t e C[z], a n ± 0). 




1. The Conjecture of Shadrin and Zvonkine 



t % :=?§t^ ••• e R. 



4 MATTHIAS ASCHENBRENNER 

Hence every element / of R can be uniquely written in the form 
/ = ^ fi f where fi G K for all i. 

i 

We call an element of R of the form at 1 , where ^ a e K , a monomial. We put 

||«|| := H + 2«i + 3i 2 + ■ ■ ■ + (n + l)i n + • • • 6 N, 
and we define a valuation v on R by setting 

v(f) := min e N for ^ / <= i?, u(0) := oo > N. 

Suppose from now on that K = Q[z] where z is a new indeterminate over Q. 
Shadrin and Zvonkine first introduce rational numbers a,d,d+k by the equation 



in the formal power series ring 

1 



l - v 



= l + V + ^ 2 H — (d = 0) 



= V + 3^ 2 + 7ip a + ■ ■ ■ {d = 1) 



1 - V 1-2^ 

1/2 1 +^4t=^ 2 + 6^ 3 + 25V 4 + --- (d = 2) 



1 - tP 1-2-0 1-3^ 



Using the numbers a,d,d+k (which turn out to be positive integers, see Lemma 1.2 
below) they then define a sequence (Lk)k>o of differential operators on R: abbre- 
viating the if-derivation of R by d n , set 

Lk = — (ln 1 ,n 1 +k 1 ■ • • 0>n r ,n r +k r t ni +ki ' " ' tn r +k r 9 ni • • ■ <9„ r (fc > 0). 

fciH hk r =k 

fei,...,fe r >0 

Note that the definition of Lk (as a if-linear map i? — > R) makes sense, since for 
every i, either 

tni+ki • ' ' Un+k T 9 ni ■ • ■ d nr (t l ) 

is zero or is a monomial which has valuation \\i\\ + ki + • • • + k r and which is 
divisible by t ni+ k 1 ■ ■ ■ t nr+ k r ; moreover, given j there are only finitely many i with 
< ||j||, and only finitely many ki,...,k r > and m,...,n r > such that 
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ini+fen • • • ,jn r +k r > 0. The first few terms of the sequence {Lk) are 

^ = E 



Li2 — «ni, 711+2^711+2 9 ni + — an 1 ,n 1 +ian 2 ,n2 + ltn 1 + ltn2 + ldn 1 dn 2 

111 ' "1,712 

-t<3 = 0.711,711+3^711+3 ^ni + — flrii, 711+1^712, "2+2 ^rii + l^n2+2^rii^n2 + 

Til 771,712 

y ^ ,7^1+2^712,712 + 1 ^711+2^712 + 1^711 ^712 ^ 
711,712 

^ ' a 7li,7li+l07l 2 ,7l 2 + l a 7l3,™3 + l ^7li +1 ^712 + 1 + 1 ! 
711,712,713 

and in general we have 

Lk = E a «i,«i+*;^i+fc <9ni + higher-order operators (fc > 0). (1.2) 

Tli 

To streamline the notation we set Lq := id/j. The argument above shows that 
for every / e R we have v(Lk(f)) ^ k + v(f), hence the sequence (z k Lk(f))k is 
summable in R. Thus one may combine the to a i^-linear map L : R — > R with 

= z kL k{f) =f + zLi{.f) + z 2 L 2 (f) + ■■■ for all / G i?. 

The operator L is used in [31] to perform a change of variables in a certain formula 
for Hurwitz numbers coming from [15]. The following proposition is established 
in [31, Proposition A. 8]. (The formula for 1^ given in [31] mistakenly omits the 
summation over n.) 

Proposition 1.1. There are rational numbers a n , n +k such that, setting 

d n (k > o) 

77 

and 

I = zh + z 2 l 2 H , 

we have L — cxp(Z), i.e., 

L w = E ^ in (/) f° r ever y f e R - 

77 

(To see that the definition of lk and I makes sense argue as for Lk and L above; 
since v(l(f)) ^ v(f) + 1 we have v(l n (f)) ^ v(f) + n for all n, hence the sum on 
the right-hand side of the equation in (1.3) exists in R.) 

After proving this proposition, Shadrin and Zvonkine make the following con- 
jecture about the form of the a n , n +k- (Again, we correct a typo in [31]: in Conjec- 
ture A. 9 replace ingjf^ by t n+k -£-.) 

Conjecture. For all k > and all n, 

(n + k + 1> 
a7i,7i+fc — c/c+i I £ ^ 

where (ck)k^i is a sequence of rational numbers, with the first terms given by (S). 
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The first step in our proof of this conjecture is to realize is that the ad,d+k are 
essentially the Stirling numbers of the second kind. We extend the definition of 
ad.d+k by setting a d d ■= 1 for every d. 



Lemma 1.2. For every d and k, 



ad,d+k 



d + k + 1 
d+1 



Proof. We expand the left-hand side of (1.1) in powers of ip: 

Ef/,)- — Vh N ' J''' I 



6=1 



d\ 



1 = WlV(-iY*-n- 1 f d 
l-bxj} ^ \d\ ; U - 1 



Now we focus on the coefficient of ip z in the last sum. By the Binomial Theorem, 
this coefficient can be written as 



' 6=0 

It is well-known that 



and 



)^»'-gG)&S ( -""(0")- 



E 

J=0 



i + l 
d+1 



a 



(1.4) 



(See, e.g., identities (6.19) respectively (6.15) in [18].) The lemma follows. 

By (1.2) and the above lemma we therefore have 

(d + k + l\ 

Lk\td) = a d ,d+ktd+k = { d+i ( td + k 

and hence 

Moreover, by definition of Ik we have lk{td) — &d,d+ktd+k for all d and fc > 0, hence 

K^d) = ^ Cid,d+kZ h t d +k 
k>0 

and thus for every n > 0: 

l n (td) = ^ Ct d ,d+kt ■ ■ ■ Cid+k 1 +---+k n - 1 ,d+k 1 + ---+k n Z kl + '" +kn td+k 1 + ---+k n - 
fei,...,fe„>0 

This yields 



exp(i)(i d ) 



E 



E ; 

fci+-+fe„=fc 

\n>0, fei,...,fc„>0 



rttd,d+fei • • • ttd+fciH hfc n _i,d+fe 



Z fe t d+ fc 



/ 
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and therefore, by (1.4) and Proposition 1.1: 



\ d+i } 



E 

feiH hfc„=fc 

n>0,fci,...,fc„>0 



1 



;a<i,d+k 1 ad+k 1 ,d+k 1 +k 2 ■ • ■ ctd+k 1 +---+k n _ 1 ,d+k- (1.5) 



It is suggestive to express this equation as an identity between matrices. We define 
{i} := ^ f° r * ^ J' ano - com bine the Stirling numbers of the second kind into a 
bi-infinite upper triangular matrix: 



/I 



S = (Sij) 





1111 



1 





1 

15 
25 
10 
1 



V 



where Si. 



(1.6) 
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We also introduce the upper triangular matrix 

/o 



A = (atj) 



V 

Then (1.5) may be written as 
Si+ij+i = ^ a 

feiH \-k n =j-i 

n>0, fei,...,fc„>0 



1 


i 


2 


n 


2 


2 


5^ 


12 


3 


-2 


2 


-4 





6 


-5 


15 
2 







10 


10 









-15 












where a,j := for i > j. 



i,i+k 1 Oti+k 1 ,i+k 1 +k 2 ' ' ' a i+fci + ---+fe„-i,j 



(i < 3) 



n=l 



or equivalently, writing 5 + := (^+1^+1)^^- and employing the matrix exponential: 

1 



n>0 



—A n = cxp(A) 



Therefore, in order to prove the conjecture from [31], we need to be able to express 
the matrix logarithm of S + in some explicit manner. We show how this can be 
done (and finish the proof of the conjecture) in Section 7 below; before that, we 
need to step back and first embark on a systematic study of a class of matrices 
(iteration matrices) which encompasses S and many other matrices of combinatorial 
significance (Sections 2 and 3), and of their matrix logarithms (Sections 4 and 5). 



2. Triangular Matrices 



In this section we let if be a commutative ring. 
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The if -algebra of triangular matrices. We construe if NxN as a if -module 
with the componentwise addition and scalar multiplication. The elements M = 
(Mij)ij e ^ of if NxN may be visualized as bi-infmitc matrices with entries in if: 



M 



/M 00 M i M 02 ••• ^ 

M 10 M n M 12 ••• 

M 20 M 2 i M 22 ••• 

V ; : •• / 



We say that M = (My) e if NxN is (upper) triangular if My = for all i, j G N 
with z > j. We usually write a triangular matrix M in the form 



/M 



M 



00 



M01 



M 02 

Mi 2 

M22 



M 03 

M13 

M23 
M33 



V 



Given triangular matrices M = (My) and M — (My), the product 

M-M:= (E fe ^M fej ) 

makes sense and is again a triangular matrix. Equipped with this operation, the 
if-submodule of if NxN consisting of all triangular matrices becomes an associative 
if-algebra ttK with unit 1 given by the identity matrix. If if is a subring of a 
commutative ring L, then tr^ is a if-subalgebra of the if-algebra ttL- We also 
define 

[M, N] := MN - NM for M, N e tt K . 

Then the if -module ttK equipped with the binary operation [ , ] is a Lie if -algebra. 

For every n we set 

tr r ^ := {M = (My) e tt K : My = for all i,j e N with i-j + n^l}. 

We call the elements of tt K strictly triangular. It is easy to verify that the 
sequence (tt K ) of if-submodules of ttK is a filtration of the if -algebra ttK, i.e., 

(1) tt° K = tt K ; 



(2) tt K D tt£- +i for all n; 

(3) teg tt£ C tt™+" for all 



and 



(4) n„tt^ = {o}. 

Clearly ttK is complete in the topology making ttK into a topological ring with 
fundamental system of neighborhoods of given by the tt K . 

The group tt K of units of tr^ has the form 

tt K = Dk x (1 + tt K ) (internal semidirect product of subgroups of tt K ) 
where Dk is the group of diagonal invertible matrices: 

D K :={M = (My) £ tt K ■ Ma e if x and My = for i ^ j}. 
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Diagonals. We say that a matrix M = (My) G txK is n-diagonal if My = 
for j ^ i + n. We simply call M diagonal if M is 0-diagonal. Given a sequence 
a = (fli)i^o S if N , we denote by diag n a the n-diagonal matrix M = (My) G if NxN 
with M j; j + „ = ai for every The sum of two n-diagonal matrices is n-diagonal. 
As for products, we have: 

Lemma 2.1. Let M = diag m a be m-diagonal and N = diag rl b be n-diagonal, 
where a — (a^), b — (bi) E K N . Then M ■ N is (m + n)-diagonal, in fact 

M ■ N = diag m+ „(a 4 • b l+m ) l>0 - 

Therefore [M,N] is (m + n) -diagonal, with 

[M, N] = diag m+n (a t • b l+m - b, ■ a l+ „) 4>0 , 

and for each k, the matrix M k is km-diagonal, with 

M k = diag fcm (ai • a i+m ■ ■ ■ a i+(fe _i) TO )^o- 

Exponential and logarithm of triangular matrices. In this subsection we 
assume that K contains Q as a subring. Then for each strictly triangular matrix 
M, the sequences (^r) n>0 an d (( — l) rl+1 ) ra>1 arc summable, and the maps 

M™ 

tr^ -> 1 + Ivk : M ^ cxp(M) := ^ — - 

n>0 U ' 



and 



1 + tx 1 K -)• tcjf : M ^ log(M) := V (-l) n+1 ^ — 112 



/ , \ 

arc mutual inverse; in particular, they are bijective. If M G tt^-, n > 0, then 
exp(M) 61 + tt£- and log(l + M) G ix n K . It is easy to see that 

cxp(M) cxp(iV) = exp(M + N) for all M, N G tt^ with MAT = 7VM. (2.1) 

In particular 

exp(M) fe = exp(kM) for all M G tt^ , fc G Z. 
We also note that given a unit C7 of ttR-, we have 

exp(C/MC/- 1 ) = C/exp(M)C/- 1 for all M G tcJr 

and 

log (UMU- 1 ) = Ulog(M)U~ 1 for all M G 1 + tr^. (2.2) 

Given M = (My)^- E ttif we define M + := (M i+ ij + i)ij E tVK- It is easy to see 
that M i ^ M+ is a if- algebra morphism tr x tr^ with M G tr^ => M+ G tr^-. 
Thus, for M e tr^ : 

exp(Af+) = exp(M)+, log(l + M+) = log(l + M)+. (2.3) 

From Lemma 2.1 we immediately obtain, for all M — diag : o where a — (aj) G if N : 

(expM)y = — — a, • a i+1 ■ ■ ■ a,j_i for all i, j G N with i < j. (2.4) 
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Derivations on the if-algebra of triangular matrices. Let 3 be a derivation 
of K, i.e., a map 3: K — > K such that 

d(a + b) = 3(a) +3(6), d(ab) = d(a)b + ad(b) for all a, b G K. 

Given M = (My) 6 ix K we let 

3(M) := (3(M„-)) S tcjf . 

Then M i-> 3(M) : it^ — > tt/f is a derivation of it^, i.e., 

3(M + AT) = 3(M) + 3(iV), d(MN) = d(M)N + Md(N) for all M, iV G tt K . 

Note that 3(tt£-) C tr^ for every n. 

We now let t be an indeterminate over K, and we work in the polynomial ring 
K* = K[t] and in the if "-algebra it if* (which contains Ixk as a K -subalgebra) . 
We equip K* with the derivation 4r. The following two elementary observations 
are used in Section 4. Until the end of this subsection we assume that K contains 
Q as a subring. 

Lemma 2.2. Let M G tr^ . Then 

4- exp(iM) = exp(tM)M. 
at 

Proof. We have (tM) n = t n M n for every n, hence 
and thus 



n! ' n! 

n>0 n>0 



- ex P (tM) = £ - ^_ J = £ _ - = ex P (^)M. 

n>0 v 7 n>0 v y 

(Similarly, of course, one also sees Jj exp(iM) = Mexp(iM), but we won't need 
this fact.) □ 

The following lemma is a familiar fact about homogeneous systems of linear 
differential equations with constant coefficients: 

Lemma 2.3. Let M, Y G tr^ and Y G tr^, . Then 
dY 

— = YM and Y\ t=Q = Y ^ Y = Y exp(tM). 

Proof. Lemma 2.2 shows that if Y = Y exp(iM) then £jL = M, and clearly 
Y\ Q = Y o exp(0) = Y . Conversely, suppose ^ = YM and Y\ Q = Y . Then 
Yi := y - Y exp(fM) G tt^, satisfies ^ = YiM and yi| t=0 = 0; hence after 
replacing Y by Yi we may assume that = YM an d ^| t _ = 0j an d need to 
show that then Y = 0. For a contradiction suppose Y ^ 0, and write Y = (Yy) 
where Yy G K* and M = (M io ) where My G if. Since Y| t _ = 0, for each i, j 
such that Yij ^Owe can write Yij — t nij Zij with G N, > 0, and Z^ G if*, 
Zij(0) 0. Choose i, j so that is minimal. Then by ^ = YM we have 

ru :i l'"> % + t»« ^ = ^ = J] Y ife M fej = 53 r»z ik M kj , 

k Y ik ^0 
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13 



dt 



and hence Zij(0) = 0, a contradiction. So Y — as desired. □ 



3. Iteration Matrices 

Let K be a commutative ring containing Q as a subring. Let A = Q[j/i, y%, . . . ] 
where (y n ) n ^i is a sequence of pairwise distinct indeterminates, let z be an inde- 
terminate distinct from each y n , and let 

Then, with x another new indeterminate, we have in the power series ring A[[x, z]]: 

where Bij = Bij(y\, y2, . ■ . ) are polynomials in y 2 , • • • ], known as the Bell 

polynomials. A general reference for properties of the is Comtet's book [11]. 
(Our notation slightly differs from the one used in [11]: B^ — Bjj.) We can obtain 
B^ by differentiating (3.1) appropriately and setting x = z = 0: 

1 d i d j , 1 d? , 

Bij — tt » ■ » explx ■ v) = —,-r-y , 

hence 

In particular, we immediately see that i?oj = and B\j = yj for j ' ^ 1. Since 

1 2' 

= S/i— r + terms of higher degree (in z) 
i) i) 

we also see that B^ — whenever i > j and -Bjj = for all j. It may also be 
shown (see [11, Section 3.3, Theorem A]) that B^ G Z[yi, . . . , and By- 

is homogeneous of degree i and isobaric of weight j. (Here each yj is assigned 
weight j.) Given a power series / G zif[[z]], written in the form 

f = ~y]f n — (/„ G if for each n ^ 1), 
n! 

n>l 
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we now define the triangular matrix 

[/] : = ([f]ij)ij eN = {Bij{h,h, • ■ • , fj-i+i)) i jeN = 

(\ 

/l /2 h h fb 

fl 3.A/ 2 4A/3 + 3/! 5^/4 + 10/2/3 
/1 6A 2 / 2 IO/2/3 + 15/i/f 

A 4 lO/f/2 
/1 



V 



g tr K 



More generally, suppose f2 = (f2„) is a reference sequence, i.e., a sequence of 
non-zero rational numbers with flo = ^i = 1. Then we define the Bell polynomials 
with respect to fi by setting 



and expanding 

ij ~ ij 

B\j = Vj f° r J ^ Ij as weu as B\j = whenever i> j and = for all j. For 
/ = XI /n finz" G zK[[z}} (/„ G K for each n ^ 1), 



(3.2) 

where B\j = (t/i, yi, . . . ) G Q[j/i, 2/2, • ■ •]• As above, one sees that B^ = and 



we define 

L/f := ([/]") i;jeN e tr* where := B« (/i,/ 2 , • • • , 

Thus, denoting the reference sequence (1/n!) by $, we have = Bij for each z, 
j and [/]* = [/] for each / e zif [[z]]. Note that by (3.2) we have, for all reference 
sequences f2, f2: 



n 



a, 



7v m - g [/] 



for all i, j, 



a L " J ' J Qi 
that is, 

(D n )~ 1 [f] n D n = (D n )- 1 [ff D n 
where D n is the diagonal matrix 



(3.3) 
(3.4) 



/fi 



V 



^2 



G tt« • 



In particular, for every reference sequence we have, with 1 denoting the constant 
sequence (1, 1,1,...): 



[ff 1 = D (D ) [/] £>*(£>") = £>" [Z] 1 (D ) . 



(3.5) 
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As first noticed by Jabotinsky [20, 21], a crucial property of [ ] n is that it converts 
composition of power series into matrix multiplication [11, Section 3.7, Theorem A]: 

[/ ° 9f = Iff ■ [gf for all f,g G zK[[z]]. (3.6) 
To see this, repeatedly use (3.2) to obtain 

]T[/ o g]% n jZ * = iU(f o gf = ilif o g = n k9 k 

j^O k^O 

= E^S ^ = E ( E^SW I 

and compare the coefficients of zK The matrix [/] n is called the iteration matrix 
of / with respect to in [11]. (To be precise, [11] uses the transpose of our [/] n .) 
For [/], the term convolution matrix of / is also in use (cf. [22]), and [/J 1 is called 
the power matrix of / in [30] . 

The subset zK x + z 2 K[[z]] of z/f[[z]] forms a group under composition (with 
identity element z), and / i-> [f] n restricts to an embedding of this group into the 
group tt£ of units of ttK- (In particular, [z] n = 1 for each f2.) As in [11], we say 
that / G zAT[[z]] is unitary if /i = 1. The set of unitary power series in K [[z]] is 
a subgroup of zK x + z 2 K[[z]} under composition, whose image under / n> [f] n is 
a subgroup of 1 + tc^ which we denote by MS^. If f2 is clear from the context, we 
simply write Mk = By (3.5), the matrix groups M^, for varying Q, are all 

conjugate to each other. We call the group of iteration matrices over K 
with respect to 0. 

Given / e AT[[z]] of the form / = z + z n+1 g with n > and g e AT[[z]] such that 
g(0) ^ 0, we say that the iterative valuation of / is n; in symbols: n = itval(/). 
(See [13].) It is easy to see that for / g Z-K]^]] and n > 0, we have / g z+z n+1 AT[[z]] 
if and only if [f] n g 1 + tr^-. For each n > we define the subgroup 

M% n := M% n (1 + teg:) = {[ff :fez + z n+1 K [[z]}} 

of M%. Then 

Mk = M n K x D M n / D • • • D M n K n D ■ ■ ■ and f| M n K n = {1}, 

n>0 

and if / g zif [[z]] is unitary with / ^ z, then n = itval(/) is the unique n > such 
that [f\ a eM% n \M% n+1 . 

As shown by Erdos and Jabotinsky [16], iteration matrices can be used to define 
"fractional" iterates of formal power series. Let t be a new indeterminate and 
K* = K[t\. 

Proposition 3.1 (Erdos and Jabotinsky). Suppose K is an integral domain, and 
let f g zAT[[z]] be unitary. Then there exists a unique power series /M g z.ff*[[z]] 
such that, writing /M := f [t] \ t=a g zif[[z]] for a g K: 

(1) /["] = z; 

(2) /[ a+1 l = /M o / /or a,beK. 
The power series is given &y 

/ [t] = E M : = E (!) ([/] - !)" e fe. . 

J>1 J ' V ' 
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Here for every n as usual ( ) = -^jt(t — 1) • ■ • (t — n + 1) 

Proof. Since [/] — 1 G tr^, the sum denning M exists in ttff*, and M\ 
for every n, by the binomial formula. Let /°* 



[/]" 

J2i>i Mij^j, and for an element a 



in a ring extension of K* write f oa := f ot . Then [/ 



Jij = M ij| t =„ = U/J'"Jij 
for every j ' ^ 1 and thus f on is the nth iterate of /: f° n = f o f o • • • o / (n times). 
In particular / 0l = / and J°( m +") = f om /°™ for all m, n. Hence if s is another 
indeterminate, then /°( s +*) = / os o / ot (in K[s, t][[z]]), since the coefficients (of 
equal powers of z) of both sides of this equation are polynomials in s and t with 
coefficients in K which agree for all integral values of (s,t). This shows that f ot 
satisfies conditions (1) and (2) (with f°' replacing everywhere). If G 
is any power series satisfying (1) and (2), then = f on is the nth iterate of /, 
for every n, and as before we deduce /W = f ot . a 

The power series (a G K) in this proposition form a subgroup of Z-ftT[[2:]] under 
composition which contains /; they may be thought of as "fractional iterates" of 
/. (This explains the choice of the term "iteration matrix.") 

Some examples of iteration matrices are collected below. Many more (in the case 
where SI = $) are given in [22]. 



([/Di 



Example. Suppose / = j^. Then 
[fh = 



1 i\ 



G N (i > 0) 



are the Lah numbers; here and below we set (^) := for i > j. (See [11, 
Section 3.3, Theorem B].) Thus if fl n = i for each n > 0, then by (3.3) 



ml 



[f]i 



hence 



/I 



Iff - 





1 2 
1 





3 4 



for i > 0, 





5 

10 
10 
5 
1 



(3.7) 



is Pascal's triangle of binomial coefficients (except for the first row). 
Example. The Stirling numbers of the second kind have the egf 



*(e*-i) 



E 

hi 



7 ■ Z J 



cf. [11, Section 1.14, (III)] or [18, (7.54)]. Hence by (3.1) we have 

[e z -l] = S, 



(3.8) 
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where S is as in (1.6). The matrix S is a unit in ttz, and it is well-known (see [11, 
Section 3.6 (II)]) that the entries of its inverse 



= (Sy 1 ) 



/l 

1-12-6 24 

1 -3 11 -50 

1 -6 35 

1 -10 
1 

V 



(3.9) 



are the signed Stirling numbers of the first kind: = (— 1) J 1 [^] , where [•?] 

denotes the number of permutations of a j-element set having i disjoint cycles. 
Thus (3.6) and (3.8) yields [log(l + z)] = S' 1 . 



4. The Lie Algebra of the Group of Iteration Matrices 

Throughout this section we let K be a commutative ring which contains Q as a 
subring. We let denote a reference sequence. We need a description of the Lie 
algebra of the matrix group M.k = M*k, generalizing the one of the Lie algebra 
of from [30]. The arguments follow [30], except that we replace the complex- 
analytic ones used there by algebraic ones. 

Definition 4.1. Let h = ^2 n h n z n 6 zif[[z]]. The infinitesimal iteration ma- 
trix of h with respect to is the triangular matrix 



/0 



(hf = ((h)?,) 















3/n ^3/> 2 
4/ii 



fii 



where = -^-ihj- i+1 . 



Note that if 17, fi are reference sequences, then 

{D n )- 1 (h) n D n = (D")- 1 ^)"!)", 

in particular 



(4.1) 
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Example 4.2. For h = ^ n h n z n G zif [[z]] we have 



(h) := (h)< 



/0 

hi 



V 



o 

ft k 



2h 



1 21 



31 h 41 h 

TT ft 3 TT " 

3!o!,_ 4! 



2/*2 112/13 

'4/ii 



where (ft)i 



hj-i + i for z > 0. 



(i-1)! 

For each n we have ft, G z n+1 K [[z]] if and only if (h) n G tr^-. We define the 
if-submodule 

TO n.» : = : ftG 



of tr^-, and we set m^- := m/; so 



m 



A' 



and P| 



m 



K 



{0}. 



n>0 



If £1 is clear from the context, we abbreviate itijf = m^- and m^- = m^'™. We set 

and we write e n if the reference sequence is clear from the context. The matrix 
e„ = e^ 2 is n-diagonal; in fact 

e„ = diag„ (j£m) S m£. 

Clearly the infinitesimal iteration matrix with respect to il of a power series from 
zK [[z]] can be uniquely written as an infinite sum 

hieo + hi&\ + • • • where h n G K for every n > 0. 

Using Lemma 2.1 one verifies easily that 

[e m , e n ] = (m - n)e m+n for all m, n. 

This implies that 

m n K = Ke n + Ke n+1 + ■■■ (n > 0) 

is an ideal of the Lie if -algebra tt^ . The main goal of this section is to show the 
following generalization of a result of Schippers [30] : 

Theorem 4.3. Let n > 0. Then exp(m^-) = M.\ (and hence log(A4^) = m^-). 

Example 4.4. Let / = j 2 -^ G zQ[[z]], and suppose 0„ = - for every n > 0. Then 
by (2.4) and (3.7) one sees easily that 



log[/] n = dia gl (0,2,3,4,...; 



/0 
2 
3 
4 


V 



/„2\fi ,- „,1 

{z } G m„. 
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We give the proof of this theorem after some preparatory results. Below we let 
f be a new indeterminate and K* = K[t\. 

Lemma 4.5. Let f e zif*[[z]] and h e zK[[z]} satisfy 

dt ~ dz 



Then 



j t [ft = [ff(hf. 



Proof. We need to show that for all i we have 

j t [f}% = ([/]»% for each j. 
For i = this is an easy computation, so suppose i > 0. We have 

and hence 

>fl h = E ( E k[fikfn k h j - k+1 ) z j . 



dz 

j^O \k=l 



Moreover 



dt ^ dt 1 llJ 3 



By the hypothesis of the lemma 



ot ot Oz Oz 



hence 

|wS = E*[/*]°^-^i = (w n w n )« 
fe=i 



3 



for each j as required. □ 

This lemma is used in the proof of the following important proposition: 
Proposition 4.6. Let h e z n+1 K [[z]], where n > 0, and set 

ft := ]T (cxp tihf^Slj z j Ez + z n+1 K*[[z}}. 

Then 

and hence 

[f t f = exp t(h) n . (4.3) 

Proof. By Lemma 2.2 we have 

^expt(h) n = (ai P t(h) n )(h) n . 
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Hence 

= J2 ((exptW )^) )^-^ 

j^l \i=l ) 

(3 \ . Q f 

^(exp t(h) n )uihj-i + i£lij z J = h. 

By Lemma 4.5 this yields f t [ft] n = [ft\ n {h) n . This shows that both Y = [f t ] n and 
Y = exp t(h) n satisfy ^ = Y (h) n and Y\ t=Q = 1. Hence [/ t ] n = cxp i(/i) n by 
Lemma 2.3. □ 

The equation (4.2) is called the formal Loewner partial differential equation in 
[30]. The following corollary, obtained by setting t = 1 in (4.3) above, shows in 
particular that exp(m^) C M.\ for each n > 0: 

Corollary 4.7. Le£ /i g [[z]], w/iere n > 0, and set 

/ := ^ (cxp (/i) n )i j O j ^ez+ * n+1 tf[[*]]. 
j>i 

2%en [f] n = exp(/i) n . 

As above we write = e| 2 . Given k\,. . . ,k n and fc = fci + • • • + k n , we have 
efei ■•■efe n =diag fe ( jy^- i(i + fci)(i + fci + k 2 ) ■ ■ • (i + fci H hfc„_i)) 

by Lemma 2.1. Now let M := (h) n where /i g Z-K"^]]. So 

M = (h) n = hieo + h 2 e 1 H 

and hence 

M n = E h ki+i- ■ ■ h k n +ie kl •••e fe „, 
fei,...,fe n 

that is, 

(M") y = 51 ft tl+ r'-/»fe„+i|j(i + fci)'--(i + fci+--'+Vi) (4-4) 

feiH yk n =j-i 

for all i, j. This observation leads to: 
Lemma 4.8. Suppose n > 0. TTien 

(M")n = (M") y = n ^._ 1) ^r% + ^i, ■ • ■ , ^-i) /or i > 2, 

where Pnj(Y , . . . , ^-2) G Q[^o ; • ■ • j ^-2] « s homogeneous of degree n and isobaric 
of weight j — 1, and independent of h. (Here each Yi is assigned weight i.) 

Proof. Set i = 1 in (4.4). Then the only terms involving hj in this sum are those 
of the form K[~ x hj ]j- j n ~ m where m 6 {!,..., n}. This yields the lemma. □ 
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An analogue of the preceding lemma (for K - 
[30]; however, the formula given there is wrong: 

Example. Suppose h — h\z + h 2 z 2 and fi = 1. Then 

/0 

/' 

M= {hy = 



and fi = l) is Lemma 3.10 of 



and hence 



M 2 = 



V 



h\ 3/iift 2 
Ah 2 



n n 


n 

u 


h 2 





2hi 2h 2 





3fti 


3/i 2 




4/ii 








2h 2 





10/ii/i 2 


6/i| 



9^? 



21/11/12 

16ft? 



V 



According to [30, Lemma 3.10] we should have, for j ^ 2: 



(M 2 



2h\hj + polynomial in hi,...,h 



However (M 2 )i2 = 3hih 2 is not of this form. 



In the proof of Theorem 4.3 we are concerned with the case where h £ z 2 JT[[;>:]], 
for which we need a refinement of Lemma 4.8: 



Lemma 4.9. Suppose h £ z 2 K[[z]) and n > 0. Then 



(M")y = 







ifn = I, 
ifl<n< j, 
otherwise. 



Proof. We have hi = 0, hence if n > 1 then (M n )ij = P$(/n, . . - ,/ij-i) by the 
previous lemma. We have M £ tr^ and hence M n £ tt^-, so (M n )y = if j — 1 < n, 
that is, if j < n. The lemma follows. □ 

Corollary 4.10. Suppose h £ z 2 K[[z}}. Then for j > 2: 



(expM)ij = — ftj 



if (&2,...,fy_l) 



where Pp(Yi, 



-,Yj-2) £ 



. ,Yj_ 2 ] is independent of h. (In particular, 



(expM)ij is polynomial in h 2 , ■ ■ ■ ,hj .) Moreover, P 2 n 
/las degree j — 1 and zs isobaric of weight j — 1. 



0, and for j > 2, P- 
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Proof. By the previous lemma we have 

(expM)y = £ ^(M n )ij = ^hj + Y ^(fti, ■ • • , Vi). 

n=l ' J n=2 

Hence 

j_1 1 

Pf(Y u . . . , ^._ 2 ) := £ -yPn^ Y l> ■ ■ ■ - ^-2) 
n=2 

has the right properties. □ 

Theorem 4.3 now follows immediately from Corollary 4.7 and the following: 

Proposition 4.11. Let f G zK[[z]] be unitary, n — itval(/). Then log[/] n G m^. 

Proof. We define a sequence {hj)j^\ recursively as follows: set h\ := 0, and as- 
suming inductively that h,2,---,hj have been defined already, where j > 0, let 

fy+i == (/j+i Pf+ifa! i- Let ft := E,->i M j e z n+1 tf[[z]] and 

M := (ft) . Then by the corollary above, we have (expM)ij = fj for every j. 
Corollary 4.7 now yields exp M = [f] n and hence log [f] n = M = (h) n G m£-. □ 

Remark. The mistake in [30, Lemma 3.10] pointed out in the example following the 
proof of Lemma 4.8 affects the statements of items 3.14 and 3.15 and the proofs of 
3.13-3.17 in loc. cit. (which concern the shape of log [/] for non-unitary / G zC[[z]]); 
however, based on the correct formula in Lemma 4.8 above, it is routine to make 
the necessary changes. For example, the corrected version of [30, Corollary 3.14] 
states that (using our notation) for ft G zC[[z]] and j^2we have 

/j . / e J hl — e hl \ 
[eMh^hj = j^j j- J +$,(fti,...,ft,-i) 

where <f>j is an entire function C 7-1 — > C. 

5. The Iterative Logarithm 

In this section we let K be an integral domain which contains Q as a subring, and 
f2 be a reference sequence. Let / G zif[[z]] be unitary. By Theorem 4.3 there exists 
a (unique) power series ft G z 2 K[[z}] such that log [f] n = {h) n . The identities (2.2), 
(3.4) and (4.1) show that ft does not depend on fi. Indeed, we have 

where ft[0] = z and h[n + 1] = h[n] of — h[n] G z" +1 i4T[[z]] for every n. 

As in [13], we call the power series ft the iterative logarithm of /, and we 

denote it by ft = itlog(/) or ft = In the following we let s, i be new distinct 
indcterminates, and we write 

/W = Y ( CX P t(f*) n )ij^J zJ e z + z" +1 if [t][[z]], n = itval(/). 

Note that does not depend on the choice of reference sequence f2. For an element 
a of a ring extension K* of K let 

/H :=/[*]| t=a £z + z" +1 ^[[z]], 
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so = z and /W = /. The notations and do not conflict with the ones 
introduced in Proposition 3.1: by (2.1) and (4.3) (in Proposition 4.6) we have 

= cxp (s + t)(h) n = cxp s(h) n • exp t(h) n = [/M] n • [/W] n - [/M o /W] n 

and hence 

yr[s+t] = f[s] Q f[t] (5 ^ 

in K[s, t][[z]]. Equation (4.2) also yields 

d f[t] 



itlog(/) 



4=0 

If a G if then = /I at l by the uniqueness statement in Proposition 3.1 and 

hence 

itlogCf^) = oit!og(/) for all a 6 K. (5.2) 
Aczel [2] and Jabotinsky [20] also showed that the iterative logarithm satisfies a 
functional equation (although [19] suggests that Frege had already been aware of 
this equation much earlier): 

Proposition 5.1 (Aczel and Jabotinsky). 



and hence 



/*•£ = /* /• ( 5 - 4 ) 

The equation (5.4) is known as Julia's equation in iteration theory. (See [24, 
Section 8.5A].) The first equation in (5.3) is simply (4.2). To show the second 

equation = /* o /M, simply differentiate (5.1) with respect to s: 



du 



u=s+t du 



d{s + t) = df[ s+t] = 9(/M /M) = dfW o 

u=s+t ds ds ds ds 

Setting s = yields the desired result. 

Suppose now that K = C. Even if / is convergent, for given a G C the formal 
power series is not necessarily convergent. In fact, by remarkable results of 
Baker [7], Ecalle [14] and Liverpool [27], there are only three possibilities: 

(1) has radius of convergence for all a G C, a ^ 0; 

(2) there is some non-zero a\ G C such that /I a l has positive radius of conver- 
gence if and only if a is an integer multiple of a\; or 

(3) has positive radius of convergence for all a G C. 

If (3) holds, then one calls / embeddable (in a continuous group of analytic 
iterates of /). This is a very rare circumstance; for example, Baker [6] and Szckercs 
[33] showed that if / is the Taylor series at of a meromorphic function on the 
whole complex plane which is regular at 0, then / is not embeddable except in the 
case where 

1 — cz 

In this case, itlog(/) = cz 2 by Example 4.4 and (5.2). Erdos and Jabotinsky [16] 
showed that in general, / is embeddable if and only if /* = itlog(/) has a positive 
radius of convergence. (See also [23, Theorem 9.15] or [29] for an exposition.) As a 
consequence, very rarely does /* have a positive radius of convergence. (However, 
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Bcalle [12] has shown that /* is always Borel summable.) In particular, we obtain 
a negative answer to the question posed in [30, Question 4.3]: if f is convergent, 
is /* convergent? Contrary to what is conjectured in [30], the converse question 
(Question 4.1 in [30]), however, is seen to have a positive answer: if /* is convergent, 
then f is convergent. 

In the next section we discuss when iterative logarithms satisfy algebraic differ- 
ential equations. 

6. Differential Transcendence of Iterative Logarithms 

Before we state the main result of this section, we introduce basic terminology 
concerning differential rings and differential polynomials. 

Differential rings. Let R be a differential ring, that is, a commutative ring R 
equipped with a derivation 3 of R. We also write y' instead of d(y) and similarly 
instead of 3™(y), where 3™ is the nth iterate of 3. The set Cr := {y £ R : y' = 0} 
is a subring of R, called the ring of constants of R. A subring of R which is closed 
under 3 is called a differential subring of R. If R is a differential subring of 
a differential ring R and y £ R, the smallest differential subring of R containing 
R U {y} is the subring R{y} :— R[y, y' , y" , . . . ] of R generated by R and all the 
derivatives y^ n ' of y. A differential field is a differential ring whose underlying ring 
happens to be a field. The ring of constants of a differential field F is a subfield of 
F. The derivation of a differential ring whose underlying ring is an integral domain 
extends uniquely to a derivation of its fraction field, and we always consider the 
derivation extended in this way. If R is a differential subring of a differential field 
F and y £ F x , then R y := {a/y n : a £ i?, n > 0} is a differential subring of F. 

Differential polynomials. Let Y be a differential indeterminate over the differ- 
ential ring R. Then R{Y} denotes the ring of differential polynomials in Y over 
R. As ring, R{Y} is just the polynomial ring R[Y, Y' , Y" , . . . ] in the distinct inde- 
terminatcs y(") over R, where as usual we write Y = Y^°\ Y' = Y^\ Y" = Y^ 2 \ 
We consider R{Y} as the differential ring whose derivation, extending the deriva- 
tion of R and also denoted by 3, is given by d(Y^) = for every n. For 
P(Y) £ R{Y} and y an element of a differential ring containing R as a differen- 
tial subring, we let P{y) be the element of that extension obtained by substituting 
y, y', . . . for Y, Y', ... in P, respectively. We call an equation of the form 

P(Y) = (where P £ R{Y}, P ^ 0) 

an algebraic differential equation (ADE) over R, and a solution of such an 
ADE is an element y of a differential ring extension of R with P(y) = 0. We say 
that an element y of a differential ring extension of R is differentially algebraic 
over R if y is the solution of an ADE over R, and if y is not differentially algebraic 
over R, then y is said to be differentially transcendental over R. Clearly to 
be algebraic over R means in particular to be differentially algebraic over R. 

Being differentially algebraic is transitive; this well-known fact follows from basic 
properties of transcendence degree of field extensions: 

Lemma 6.1. Let F be a differential field and let R be a differential subring of F. 
If f £ F is differentially algebraic over R and g £ F is differentially algebraic over 
R{f}, then g is differentially algebraic over R. 
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Differential transcendence of iterative logarithms. Let now K be an inte- 
gral domain containing Q as a subring, and let z be an indeterminate over K . We 
view K l\z]} as a differential ring with the derivation The ring of constants of 
K [[z]] is K. We simply say that / G -K"[[z]] is differentially algebraic or differ- 
entially transcendental if / is differentially algebraic respectively differentially 
transcendental over K[z\. If / G ^[[z]] is differentially algebraic, then / is actually 
differentially algebraic over K, by Lemma 6.1. 

As above, we let f be a new indeterminate over K, and K* = K[t\. The goal of 
this section is to show: 

Theorem 6.2. Let f G zif[[z]] be unitary. Then /» G z 2 -K"[[z]] is differentially 
algebraic if and only if G is differentially algebraic, if and only if f^ 

is differentially algebraic over K* . 

Before we give the proof, we introduce some more terminology concerning dif- 
ferential polynomials, and we make a few observations about how the derivation 
■j^ of ii"[[z]] and composition in if[[z]] interact with each other, in particular in 
connection with solutions of Julia's equation. 

More terminology about differential polynomials. Let R be a differential 
ring and P G R{Y}. The smallest r G N such that P G R[Y, Y',..., Y^} is called 
the order of the differential polynomial P. Given a non-zero P G R{Y} we define 
its rank to be the pair (r, d) G N 2 where r — ordcr(P) and d is the degree of P in 
the indeterminate Y^ r \ In this context we order N 2 lexicographically. 

For any (r + l)-tuplc i = (i , ■ . ■ , i r ) of natural numbers and Q G R{Y}, put 

Q i := Q i °{Q') il ■ ■ ■ (Q^) ir . 

In particular, Y i = Y ia {Y') il ■ ■ ■ (Y^) ir , and y i = y io (y') 11 ■ ■ ■ (y {r) ) lr for y G R. 

Let P G R{Y} have order r, and let i — (i , ■ ■ ■ , i r ) range over N 1+r . We denote 
by Pi G R the coefficient of Y % in P; then 

P{Y) = Y,PiY i - 

i 

We also define the support of P as 

suppP := {i : Pi ^ 0}. 

We set 

\i\ := io H Mn ||*|| : = h + 2 *2 H Yri r . 

For non-zero P G R{Y} we call 

dcg(P) = max \i\, wt(P) = max ||z|| 

iGsuppP iGsuppP 

the degree of P respectively weight of P. We say that P is homogeneous if 

\i\ = dcg(P) for every i G suppP and isobaric if ||i|| = wt(P) for every i G suppP. 

Transformation formulas. Let X be a differential indeterminate over -K"[[z]]. An 
easy induction on n shows that for each n > there arc differential polynomials 
G mn G Z{X} (1 ^ to < n) such that for all / G and h G we have 

(h (n) o /) • (Z') 2 "- 1 = Gm(/) (h o /)' + G 2n (/) (ft o /)" + ■ • ■ + G„„(/) (ft o /)<"). 



24 



MATTHIAS ASCHENBRENNER 



Moreover, G m n has order n — m + 1, and is homogeneous of degree n — 1 and 
isobaric of weight 2n — m — 1. Set G mn :=Oifm>iiorm = 0<n, and 
Goo : = e Then the G mn satisfy the recurrence relation 

= (1 - 2n)G mn X" + (G' mn + G m -x, n )X' (m > 0). 

Organizing the G mn into a triangular matrix we obtain: 

({x')- 1 

1 -X" 3(X") 2 - X'X^ 
X' -3X'X" 

(X') 2 



G := (G 



mnjm.n 



V 



(6.1) 



Note that G n 



(X> 



for every n. Now set 



H kn = (™)x (m - k+1) Gmn G Z{X} 



for fc = 0, 



,n. 



So if we define the triangular matrix 



B := (B, 



km) 



/X' X" x& x& 
X' 2X" 

X' 3X" 
X' 



V 



where B km = ^j X (m - k+1 '> for m > fc, 



then 



# ; — {Hkn) — B ■ G — 

(I X" X'X^-{X"Y (X'fX^ -\X'X"X^ +3(X") 3 -..\ 
X'X" -3X'(X") 2 + 2(X') 2 X^ 

{X') 2 

{x'f 



X' 



V 



Each differential polynomial H kn has order at most n — k + 1, and if non-zero, is 
homogeneous of degree n and isobaric of weight 2n — fc. Note that for n > 0, -Hon 
has the form 

n 

Hon = ^ X {m+1 ^G mn = (X')"- 1 *^ 1 ) + ff„ where ff„ e Z[A", . . .,X™]; 

m—l 

in particular order(iJ n) = n + 1 > order(iJ / t„) for fc = 1, . . . , n. 
Let now / e zi-T^z]] and /i 6 -K"[[z]] satisfy Julia's equation 

h-f' = ho f. 

We assume / ^ (and hence /' ^ 0). Then for every n: 
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Let R := K{X}x>, and denote the i?-algebra automorphism of R{Y} with 
^ (X') 1 - 2 " (Hon Y + H ln ¥' + ■■■ + H nn for every n 

also by if. Then for every P £ K{Y} we have 

P(h)of = H{P)\ x=fY=h . 
Note that for every i £ N and n we can write 

(x') (2 " _1)i • h((y^y) = (x'y^-VY^x^+vy + a, 

where a t £ Z[X', . . . , X^ n+1 \ Y, Y' , . . . , with deg x< „ +1) 0i < i. 
Hence given i = (i , . . . ,i r ) £ W +1 , setting d = \i\ and w = \\i\\, we may write 

(X') 2w -< 1 ■ H(Y i ) = (X r ) w - d (X'yY d + at 

where a t £ Z[X', . . . , X {r+1 \ Y, Y' Y {r) ] with deg x(r+ i) a< < i r . 

Proof of Theorem 6.2. Let / G zK [[z]] be unitary. Suppose hrst that is 
differentially algebraic over K* . Let P £ K*{Y} be non-zero of lowest rank such 
that P(f^) — 0. Differentiating with respect to t on both sides of this equation 
yields 

i=0 

Here r = order(P) and P*(Y) £ K*{Y} is the differential polynomial obtained by 
applying ^ to each coefficient of the differential polynomial P. Now by Proposi- 
tion 5.1 we further have 

(i) 

\ / / 

Since d y P (r ) has lower rank than P, by choice of P we have gf^Tj (/'*') ^ 0. Hence /* 
satisfies a non-trivial (inhomogeneous) linear differential equation with coefficients 
from K *{/'*'}, and so by Lemma 6.1, is differentially algebraic over K*. Specializing 
t to a suitable rational number in an ADE over K* satisfied by /* shows that then 
/* also satisfies an ADE over K, that is, /* is differentially algebraic over K. 

Conversely, suppose that /» is differentially algebraic. Let P £ K {Y} be non- 
zero, of some order r, such that P(/*) = 0. Then 

H{P)(f [t \U) = P(U)of^ =0. 

Let d = deg y (r) P. By the remarks in the previous subsection, for sufficiently large 
JVeNwe have 

{X') N H{P)= Pi(X') N -^Y^+A 

i:i r —d 

where A £ K[X', . . . , X {r+1 \ Y, Y' , . . . , Y {r) ] with dcg x(r+ i) A < d. 
For such N, the differential polynomial 

Q(X) := (X') N H(P)\ Y=h £R{X} 

is non-zero, where R = K{f*}, and satisfies Q(f [t] ) = 0. Thus /M is differentially 
algebraic over R and hence (by Lemma 6.1) over K, as required. □ 
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Let J 7 be a family of elements of -K"[[z]] • Following [10] we say that T is coherent 
if there is a non-zero differential polynomial P G i4T[z]{y} such that P(f) = for 
every /e J. If J 7 is coherent, then P with these properties may actually be chosen 
to have coefficients in K; see [10, Lemma 2.1]. If J- is not coherent, then we say 
that T is incoherent; we also say that T is totally incoherent if every infinite 
subset of T is incoherent. From the previous theorem we immediately obtain a 
result stated without proof in [10]: 

Corollary 6.3 (Boshernitzan and Rubel [10]). Let f G zif [[z]] be unitary and let 
J~ '■= {.f^°\ f^K f^ 2 \ ■ ■ ■} be the family of iterates of f . Then exactly one of the 
following holds: 

(1) /* is differentially algebraic and T is coherent; 

(2) /» is differentially transcendental and F is totally incoherent. 

Proof. By the theorem above, it suffices to show: if is differentially algebraic, 
then T is coherent, and if /M is differentially transcendental, then T is totally 
incoherent. The first implication is obvious (specialize t to n in a given ADE 
for /'*]). For the second implication, suppose T is not totally incoherent. Then 
there exists an infinite sequence (nA of pairwise distinct natural numbers such that 
{/["•]} is coherent. Let P G K{Y}, P ^ 0, be such that P(f [n ' ] ) = for every L 
With g := P(f^) G if* [[2]] we then have g\ = for every i; thus g — (since 
the coefficients of g are polynomials in t with coefficients from the integral domain 
K* of characteristic 0). This shows that is differentially algebraic. □ 

7. The Iterative Logarithm of e z - 1 

In this section we apply the results obtained in Sections 4 and 5 to the unitary 
power series / = e z — 1 G zQ[[z]]. Recall that the iteration matrix [e z — 1] of this 
power series is the matrix S — (Sij) G 1 + itq consisting of the Stirling numbers 
Sij = {{} of the second kind (cf. (3.8)). 

Proof of the conjecture. We first finish the proof of the conjecture stated in 
Section 1. The matrix S is related to A = (ay) G itq via the equation 

S+ =exp(A), 

or equivalently (cf. (2.3)): 

A = Iog(5)+. 

(Recall: for a given matrix M = (My) G ttQ we defined M + = (Mi+ij+i)ij G ttq.) 
The conjecture postulates the existence of a sequence (c„)„j>i of rational numbers 
such that 

a ij =c j -i +1 [ . ) fori<j. (7.1) 



1 

This now follows easily from the results of Section 4: 

Proposition 7.1. Let h = itlog(e z — 1) G z 2 Q[[z]] ; write h = J2n^i^nZ n where 
h n G Q, and define c„ := n\ h n for n ^ 1. Then (7.1) holds, and 

E(-l) fe+1 \n 2 \ (n 3 \ f n k 

l<ni<-<»ii_i<TH=n 

for every n > 1. 



k ImJ [n 2 J tn fc _! 
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Proof. We have log (5) = (h) by Theorem 4.3. Hence, using the formula for (h)ij 
from Example 4.2 we obtain for i < j, as required: 



(j + 1)1 

o-ij = (h}i+ij + i = —r hj-i+i 



+ 



i\ J ^ i\(j-i + iy. 

The displayed identity for c„ follows from c„ = {h)\ n = log(S)i„. 



3 + 1 

i 



□ 



We note that the c„ may also be expressed using the Stirling numbers of the 
first kind, using (h) = — \og(S" 1 ): 



E 



(-i) 



k+n—m 



Kni<---<n k - 1 <n k =n 



n k -l 



(n> 1). 



Proof of the convolution identity. We now turn to the convolution identity (C) 
for Stirling numbers stated in the introduction. Jabotinsky's functional equation 
(5.4) for / = e z — 1, writing again h — /*, reads as follows: 

h o (e z - 1) = e z h. 

Taking derivatives on both sides of this equation and dividing by e z we obtain: 



h'o(e z -l) = h + ti. 



(7.2) 



Now define, for M E 1 + ttk: 



C_1 )n 

A(M) := V [ — 1 T (M - l) n £ 1 + t^, 
' n + 1 



so 



A(M) ■ (M - 1) = log(M). 
For later use we note that then for every j > 1 : 



(7.3) 



5^A(M)i fc M fcij+ i = 5^A(M)i fc (M-l) fciJ - + i = 
fe=i fe=i 

(A(M) • (M - = log(M) M+1 , (7.4) 

where in the last equation we used (7.3). 
Taking M — S we compute 



log(S) 



1 -\ 
3 












1 


2 


11 


2 


5 5 


12 


-2 


2 


-4 


6 


-5 


15 
2 





10 


10 







-15 
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and 



/l 
1 



A(5) 



_ i 

2 

1 - 












1 


2 


n 


2 


5 3 


12 


3 




25 


2 


2 




1 


-3 


15 6 
2 




1 


-5 

1 



V 



We observe that the first row of A(S) agrees with the first row of log(S) shifted by 
one place to the left. (This is simply a reformulation of the formula (C).) 

Proposition 7.2. For every j > 1, 

A(S)y=]og(S) llj+1 . 

Proof. As observed in (7.4), 



k=l 



.7 + 1 

k 



Cj+i for j > 1. 



(7.5) 



On the other hand, by (7.2) we have [h 1 ] ■ S = [h + h']; thus 



f 7 + ll 

E Cfe+1 { fc } = Xj [ ft V s Wi = (M • S )l,J+l = [ ft + ''li J+i = Cj'+i + c i+2 

„ j j + 1 



and hence 



k=i 



c J+ i for j ^ 1. 



(7.6) 



An easy induction on j using (7.5) and (7.6) now yields A(S)\j = cj + \ = log(5)i J+ i 
for each j ^ 1, as claimed. □ 

Differential transcendence of the egf of (c„). It is easy to see that for n > 0, 
the nth iterate <j>^ of 4> = e z — 1 is a solution of an ADE over Q of order n. 
However, it is well-known that <j>^ does not satisfy an ADE over <C[z] of order < n. 
(See, e.g., [5, Corollary 3.7].) The egf of the sequence (c„) is itlog^, hence from 
Corollary 6.3 we obtain the fact (mentioned in the introduction) that this egf is 
differentially transcendental. In fact, Bcrgweiler [9] showed the more general result 
that if / is (the Taylor scries at of) any transcendental entire function, then 
itlog(/) is differentially transcendental (equivalently, by Corollary 6.3, the family 
of iterates of / is totally incoherent). Moreover, by the results quoted at the end 
of the previous section, itlog^ is not convergent. (This can also be shown directly; 
cf. [25].) See [3] for a proof of a common generalization of these two facts. 
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